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Several physical realizations of quantum bits have been proposed, including trapped ions, NMR 
systems and spins in nano-structures, quantum optical systems, and nano-electronic devices. The 
latter appear most suitable for large-scale integration and potential applications. We suggest to use 
low-capacitance Josephson junctions, exploiting the coherence of tunneling in the superconducting 
state combined with the possibility to control individual charges by Coulomb blockade effects. These 
systems constitute quantum bits (qubits), with logical states differing by one Cooper-pair charge. 
Single- and two-bit operations can be performed by applying a sequence of gate voltages. The phase 
coherence time is sufficiently long to allow a series of these steps. The ease and precision of these 
manipulations depends on the specific design. Here we concentrate on a circuit which is most easily 
fabricated in an experiment. 

In addition to the manipulation of qubits the resulting quantum state has to be read out. This 
can be accomplished by coupling a single-electron transistor capacitively to the qubit. To describe 
this quantum measurement process we study the time evolution of the density matrix of the coupled 
system. Only when a transport voltage is turned on, the transistor destroys the phase coherence 
of the qubit; in this case within a short time. The measurement is accomplished after a longer 
time, when the signal resolves the different quantum states. At still longer times the measurement 
process itself destroys the information about the initial state. We present a suitable set of system 
parameters, which can be realized by present-day technology. 



I. INTRODUCTION 

The investigation of nano-scale electronic devices, such 
as low-capacitance tunnel junctions or quantum dot 
systems, has always been motivated by the perspec- 
tive of future applications. By now several have been 
demonstrated, e.g. the use of single-electron transis- 
tors (SET) as ultra-sensitive electro-meters and single- 
electron pumps. From the beginning it also appeared at- 
tractive to use these systems for digital operations needed 
in classical computation (see Ref. jjj for a review). Ob- 
viously single-electron devices would constitute the ulti- 
mate electronic memory. Unfortunately, their extreme 
sensitivity makes them also very susceptible to fluctu- 
ations and random background charges. Due to these 
problems — and the continuing progress of conventional 
techniques — the future of SET devices in classical digi- 
tal applications remains uncertain. 

The situation is different when we turn to elements for 
quantum computers. They could perform certain calcula- 
tions which no classical computer could do in acceptable 
times by exploiting the quantum mechanical coherent 
evolution of superpositions of states Here conven- 
tional systems provide no alternative. In this context, 
ions in a trap, manipulated by laser irradiation are the 
best studied system [|]||]. However, alternatives need to 
be explored, in particular those which are more easily em- 
bedded in an electronic circuit. From this point of view 
nano-electronic devices appear particularly promising. 

The simplest choice, normal-metal single-electron de- 
vices are ruled out, since — due to the large number of 
electron states involved — different, sequential tunneling 



processes are incoherent. Ultra-small quantum dots with 
discrete levels, and in particular, spin degrees of freedom 
embedded in nano-scale structured materials are candi- 
dates. They can be manipulated by tuning potentials and 
barriers ||. However, these systems are difficult to fabri- 
cate in a controlled way. More attractive appear systems 
of Josephson contacts, where the coherence of the super- 
conducting state can be exploited and the technology is 
quite advanced. Macroscopic quantum effects associated 
with the flux in a SQUID have been demonstrated [[|. 
Quantum extension of elements based on single flux logic 
have been suggested , and efforts are made to observe 
one of the elementary processes, the coherent oscillation 
of the flux between degenerate states [[)) . 

We have suggested to use low-capacitance Joseph- 
son junctions, where Cooper pairs tunnel coherently 
while Coulomb blockade effects allow the control of in- 
dividual charges. They provide physical realizations of 
quantum bits (qubits) with logical states differing by the 
number of Cooper pair charges on an island. These junc- 
tions can be fabricated by present-day technology. The 
coherent tunneling of Cooper pairs and the related prop- 
erties of quantum mechanical superpositions of charge 
states have been discussed and demonstrated in experi- 
ments |fl]-(l6|]. In particular, in a recent experiment |l7| ] 
Nakamura et al. have observed time-resolved coherent 
oscillations of charge in such Josephson junction setup. 

We concentrate here on the simplest design of the 
qubit, which is sufficient to introduce the ideas, and it 
is most easy for fabrication. It should be added that 
other ideas in this direction have been discussed |18| , ^9[ . 
In Ref. [jl9| an improved design has been proposed that 
relaxes requirements to the parameters of the circuit and 
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makes the manipulation procedures more flexible. In Sec- 
tion H we will introduce the system and show how single- 
and two-bit operations (gates) can be performed by the 
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application of sequences of gate voltages. In Section 
we analyze the influence of dissipation and fluctuations 
and conclude that for a proper choice of parameters the 
phase coherence time is sufficiently long to allow a large 
number of gate operations. 

In addition to the controlled manipulations of qubits, 
quantum computation requires a quantum measurement 
process to read out the final state. The requirements 
for both steps appear to contradict each other. Dur- 
ing manipulations the dephasing should be minimized, 
whereas a measurement should dephase the state of the 
qubit as fast as possible. The option to couple the mea- 
suring device to the qubit only when needed is hard to 
achieve in nano-scale systems. The alternative described 
in Section IV , is to couple a normal-state single-electron 
transistor capacitively to a qubit pcj ]. During the ma- 
nipulations the transport voltage of the SET is turned 
off, and the SET only acts as an extra capacitor. To per- 
form the measurement the transport voltage is turned 
on. The dissipative current through the transistor de- 
phases the qubit and provides the read-out signal for the 
quantum state. We describe this quantum measurement 
process by considering explicitly the time-evolution of the 
density matrix of the coupled system. We find that the 
process is characterized by three time scales: a short de- 
phasing time, the longer 'measurement time' when the 
signal resolves the different quantum states, and finally 
the mixing time after which the measurement process it- 
self destroys the information about the initial state. Sim- 
ilar nonequilibrium dephasing processes |2Uj23] have re- 
cently been demonstrated experimentally~p4f . 

In Section we discuss system parameters and sug- 
gest suitable sets which can be realized by present-day 
technology. We further compare with related work. 



II. JOSEPHSON JUNCTION QUBITS 



A. Qubits and single-bit gates 

The simplest Josephson junction qubit is shown in 
Fig. 0a. It consists of a small superconducting island 
connected by a tunnel junction, with capacitance Cj and 
Josephson coupling energy E j , to a superconducting elec- 
trode. An ideal voltage source, V^b, is connected to the 
system via a gate capacitor C (fluctuation effects will be 
discussed later). We choose a material such that the su- 
perconducting energy gap A is the largest energy in the 
problem, larger even than the single-electron charging 
energy to be discussed below. In this case quasi-particle 
tunneling is suppressed at low temperatures, and a situ- 
ation can be achieved where no quasiparticle excitations 
exist on the island*. 



E, 



T T T 



a) 



b) 



FIG. 1. Ideal a) one qubit and b) multi-qubit systems. 

In the following we will consider the situation where 
only Cooper pairs tunnel in the superconducting junc- 
tion. This system is described by the Hamiltonian 



H = 4E qh (n- 



^qbj 



Ej cos . 



(1) 



Here, n is the number (operator) of extra Cooper pair 
charges on the island (relative to some neutral refer- 
ence state) and the phase variable is its conjugate 



* Under suitable conditions the superconducting state is to- 
tally paired, i.e. the number of electrons on the island is even, 
since an extra quasi-particle (odd number of electrons) costs 
the extra energy A. This 'parity effect' has been established in 
experiments below a crossover temperature T* ~ A/lniV c ff, 
where N c g is the number of electrons in the system near the 
Fermi energy [ p^|j2^ -|2^] . For a small enough island, T* is 
typically one order of magnitude smaller than the supercon- 
ducting transition temperature. For the case that — e.g. due 
to the initial preparation — an unpaired excitation exists on 
the island, a channel should be provided for the quasiparticle 
to escape to normal parts of the system M . 
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n = —ihd/d(ti@). The charging energy of the super- 
conducting island is characterized by the scale E q i, = 
e 2 /2(C + Cj), while the dimensionless gate charge, n q b = 
CV q b/2e, acts as a control field. We consider systems 
where -E q b is much larger than the Josephson coupling 
energy, E q ^ 3> Ej. In this regime a convenient basis is 
formed by the charge states, parameterized by the num- 
ber of Cooper pairs n on the island. In this basis the 
Hamiltonian (]l|) reads 

H = J2 / LE qh (n-n qh ) 2 \n)(n\ 

n 

-l Ej (\n)(n + l\ + \n + l)(n\) . (2) 

For most values of n q b the energy levels are dominated 
by the charging part of the Hamiltonian. However, when 
n q b is approximately half-integer and the charging ener- 
gies of two adjacent states are close to each other, the 
Josephson tunneling mixes them strongly (see Fig. 0). 




1/2 n qb 

FIG. 2. The charging energy of the superconducting elec- 
tron box is shown (solid lines) as a function of the gate charge 
n q b for different numbers of extra Cooper pairs n on the is- 
land. Near degeneracy points the weaker Josephson coupling 
energy mixes the charge states and modifies the energy of the 
eigenstates (dotted line) . In this regime the system effectively 
reduces to a two-state quantum system. 

We concentrate on the voltage interval near a degen- 
eracy point of two charge states, say n = and n = 1. 
For the parameters chosen, further charge states can be 
ignored, and the system (Q) reduces to a two-state model, 
with Hamiltonian which can be written in spin-^ nota- 
tion in terms of Pauli spin matrices a — a x ,a y ,a z as 

H = -J5 c h(n q b)CT« - -Eja x . (3) 

The charge states n = and n = 1 correspond to the spin 
basis states |j) = (9) and |f) = (J), respectively. The 
charging energy -E c h(«qb) = 4i? q b(l — 2n q b), equivalent 
to a magnetic field in z-direction in the spin problem, is 
controlled by the gate voltage. For convenience we can 
further rewrite the Hamiltonian as 

H = ItAE(t]) (cost? a z - sin?7 a x ) , (4) 



where the mixing angle r\ = tan -1 [-Ej/-E c h(n q b)] deter- 
mines the direction of the effective magnetic field in the 
xz-plane, and the energy splitting between the eigen- 
states is AE(-q) = v^ch^qb) + E j = Ej/sint]. At the 
degeneracy point, r/ — tt/2, it reduces to Ej. The eigen- 
states are denoted in the following as |0) and |1). For 
some chosen value of n q b they are 

|0)= cos 1 1|) + sin \ |t) 
|l)=-sin|||) + co S ||T). (5) 

For later convenience we can rewrite the Hamiltonian 
in the basis of eigenstates. To avoid confusion we intro- 
duce a second set of Pauli matrices, p, which operate in 
the basis |0) and |1), while reserving a for the basis of 
charge states | j ) and | f ) . By definition the Hamiltonian 
then becomes H = ^AE(i])p z . 

By changing the gate voltage we can perform the re- 
quired one-bit operations (gates). If, for example, one 
chooses the idle state far away from degeneracy, the 
eigenstates |0) and |1) are close to |!) and |f), respec- 
tively. For definiteness, we choose the eigenstates |0) and 
|1) at the idle point as logical basis states. Then switch- 
ing the system suddenly to the degeneracy point for a 
time At and then switching back produces a rotation in 
spin space, 

/ a , \ /. \ /cos a isina\ , , 

L/x_ bit (Ai) = exp^) = . s . nQ cQsa j , (6) 

where a = EjAt/2h. Depending on the value of At, a 
spin flip can be produced, or, starting from |0), a su- 
perposition of states with any chosen weights can be 
reached. This is exactly the way the experiments of Naka- 
mura et al. @ were performed. At the same time, keep- 
ing the system at the idle point we achieve a phase shift 
between two logical states by the angle Ai?(ryidi e )At. A 
different phase shift during the same time period At can 
be achieved using a temporary change of n q b by a small 
amount which changes the energy difference between the 
ground and excited states. 

The example presented above provides an approximate 
spin flip for the situation where the idle point is far from 
degeneracy and E q \, ^> Ej. But a spin flip in the logi- 
cal basis can also be performed exactly. It requires that 
we switch from the idle point rjidie to the point where 
the effective magnetic field is orthogonal to the idle one, 
V = r /idic + 7r/ 2. This changes the Hamiltonian from 
H = \AE{rr, A x e )p z to H = iA£(77 idlo + tt/2) Px . To 
achieve that, the dimensionless gate charge n q b should 
be increased by i?j/(4i? q bsin27yidic). In the limit dis- 
cussed above, yyidic 1, the operating point is close to 
the degeneracy, r/ = n/2. 

Unitary p x - and p 2 -rotations described above are suffi- 
cient for all manipulations with a single qubit. By using 
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a sequence (of no more than three) such elementary ro- 
tations we can achieve any unitary transformation of the 
qubit's state. 

In our discussion so far elementary rotations should 
have been performed one immediately after another. 
However, sometimes the quantum state should be kept 
intact for a certain time interval, for instance, while other 
qubits are manipulated during the computation. Even in 
the idle state r\ = ryidic, the energies of the two eigenstates 
are different. Hence their phases evolve relative to each 
other, which leads to the quantum mechanical 'coherent 
oscillations' of a system which is in a superposition of 
eigenstates. We have to keep track of this time depen- 
dence as is demonstrated by the following example (for 
simplicity again on the approximate level). Imagine we 
try to perform a rotation by angle a in the spin space, 
which we can do by a suitable choice of At in (||) , and af- 
ter some time delay r we perform the reverse operation. 
The unitary transformation for this combined process is 

/ • \ f.AE(T] mc )r \ f. \ 
exp ( - iot(j x ) ■ exp I % — a z I • exp (iaa x ) . 

Clearly the result depends on the intermediate time r 
and differs from unity or a simple phase shift, unless 
AE(rjm c )T/2h = nix. The time-dependent phase factors, 
arising from the energy difference in the idle state, can be 
removed from the eigenstates if all the calculations are 
performed in the interaction representation, with zero- 
order Hamiltonian being the one at the idle point. In 
this way the information which is contained in the am- 
plitudes of the qubit's states is preserved. However, there 
is a price for this simplification, namely the transforma- 
tion to the interaction representation introduces addi- 
tional time dependence in the Hamiltonian during the 
operations. Thus a general 1-bit operation induced by 
switching at to from ?7;dic to r\ for some time At is de- 
scribed by the unitary transformation 

U(t , At, n) = e £^dio)(to+At) e -%H( v )At g-jj.fffiftdujto . ( 7 ) 

This demonstrates that the effect of the operations de- 
pends not only on their time span but also on the moment 
when they start. 

The transformation (Q) is the elementary operation in 
the interaction picture in the sense that it can be per- 
formed in one step (one voltage switching) . Can we per- 



form an arbitrary single-bit gate using such an opera- 
tion? To answer this question we note that all the uni- 
tary transformations of the two-state system form the 
group SU(2). The group is three-dimensional, hence one 
needs in general three controllable parameters to obtain 
a given single-bit operation. As is obvious from (0), IA 
depends exactly on three parameters: 77, to and At. It 
can be shown that any single-bit gate can be performed 
in one step if the voltage (i.e. 77) and the times to, At 
are chosen properly.^ (The dependence on to is periodic 
with the period hj AE^rjidie)- Hence the waiting time to 
a new operation is restricted to this period.) 

As a final remark in this subsection, we notice that 
our choice of the logical basis is by no means unique. As 
follows from the preceding discussion, we can perform 
x- and z-rotations in this basis, which provides sufficient 
tools for any unitary operation. On the other hand, since 
we can perform any unitary transformation, we can as 
well perform x- and z-rotations in any basis of our two- 
dimensional system. Therefore, any basis can serve as 
the logical one. The Hamiltonian at the idle point is 
diagonal in the basis (|B|), while the controlled part of 
the Hamiltonian, the charging energy favors the charge 
basis. For this reason manipulation procedures for x- 
and z-rotations, in the interaction representation, are of 
comparable simplicity in all bases. The preparation pro- 
cedure (thermal relaxation at the idle point) is easier de- 
scribed in the eigenbasis, while coupling to the meter (see 
Section [Tv| ) is diagonal in the charge basis. So, there is 
no obvious choice for the logical states and it is a matter 
of convention. 



B. Many-qubit system 

For quantum computation a register, consisting of a 
(large) number of qubits is needed, and pairs of qubits 
have to be coupled in a controlled way. Such two-bit 
operations (gates) are necessary, for instance, to create 
entangled states. For this purpose we couple all qubits by 
one mutual inductor as shown in Fig. ^b. One can easily 
see that for L = the system reduces to a series of uncou- 
pled qubits, while for L = 00 they are coupled strongly. 
Finite values of L introduce some retardation in the in- 
teraction. The Hamiltonian of this system, consisting of 



Tt is customary to consider rotations about two axes, say, 
x- and z-rotations, as device-independent elementary single- 
qubit gates. Then, any single-bit gate in a quantum algorithm 
is described as a combination of three consecutive elemen- 
tary rotations given by their angles (e.g. Euler angles a, f3, 
7). To realize such a gate O we need to solve the equation 
U(j], to, At) = 0(a,/3, 7) to find the triple of rj, to and At. 
Since usually only a few different single-bit gates are used, 
corresponding triples can be calculated in advance. 



4 



N qubits and an oscillator formed by the inductance and 
the total capacitance of all qubits is |1(J 



H 



J2iLi |4-E q b(nj - n qb ,i) 2 



2NC t 



Here $ is the flux in the mutual inductor, $o = h/2e 
is the flux quantum, and q — —iTid/d$>, the variable 
conjugated to is related to the total charge on the 
gate capacitors of the qubits. For the oscillations in 
this LC-circuit the junction and gate capacitor of each 
qubit act in scries, hence the relevant capacitance is 
C7 1 = Cj" 1 + C -1 . The voltage oscillations of the LC- 
circuit affect all the qubits equally, thus $ is coupled to 
each of the phases O;. The reduction factor Ct/Cj de- 
scribes the screening of these voltage oscillations by the 
gate capacitors. Here we have chosen to express the cou- 
pling via a shift of the phase in the Josephson coupling 
terms arising from the voltage oscillations. This form is 
reminiscent of the usual description of SQUIDs, except 
that $ is a dynamic variable and its effect is reduced by 
the ratio of capacitances. Alternatively, we could have 
added the oscillating charges to the gate charges in the 
charging energy. Both forms are equivalent and related 
by a canonical transformation. This and a more detailed 
derivation of (g) has been presented in Ref. fjcfl . 

The oscillator described by the charge q and flux $ 
with characteristic frequency — (NLCt)^ 1 ^ 2 pro- 

duces an effective coupling between the qubits. Wc 
choose parameters such that 



AE( V ) « hJ L N J 



(9) 
(10) 



The first condition (g) assures that the oscillator remains 
in its ground state at all relevant operation frequencies. 
I.e. the logical operations on the qubits are not affected 
by excited states. The second condition ( |i~f]| ) prevents the 
Josephson coupling in (]|) from being "washed out" by 
the fluctuations of <f>. Since they are limited by L, hence 
($ 2 )/L w > this condition imposes only weak con- 

straints on the parameters. 

Although the LC-oscillator remains in its ground state, 
it provides an effective coupling between the qubits. To 
analyze this we expand in ( |8|) the Josephson coupling 
terms in $ and, because of (Rw, neglect powers higher 
than linear. The linear term is /<£>, where the current 
through the inductor is given by the sum of qubits' con- 
tributions, 



C t 2itEj ^ . 
i = 75 — — ) sin &i . 



(11) 



The linear term and the unperturbed Hamiltonian of the 
oscillator can be combined to a square, 



2NC t 



$ 2 

2L 



— +7$ = 



2NC t 



(<5> + LI)'< 
2L 



As long as the frequency of the oscillator is large com- 
pared to characteristic frequencies of the qubit 's motion 
(0), one can use the adiabatic approximation and treat 
the slow qubit's variables (Oi) as constant to find the en- 
ergy levels of the oscillator. Since the lowest level of the 
first two terms in the rhs of Eq.([l2|), equal to Hujlc/2, 
does not depend on /, the correction to the ground state 
energy is given by the last term. This term provides the 
effective coupling between the qubits 



-Hint — —Ej 



where the energy scale is 



sin 0, 



(13) 



E L = 2n 



2 C t 2 EjL 



(14) 



In the spin-i notation sinOi = \<Ty and the interaction 
term becomes (up to constant terms) 



E L 



E 



£T« ff Ci) 



(15) 



The ideal system would be one where the coupling be- 
tween different qubits could be switched on and off, leav- 
ing the qubits uncoupled in the idle state and during 1-bit 
operations. This option requires a more complicated de- 
sign |is|| . With the present, simplest model the qubits are 
coupled at all times. But even in this case we can con- 
trol the coupling in an approximate way by tuning the 
energies of the selected qubits in and out of resonance. If 
El is smaller or comparable to Ej and the gates voltages 
Vqh^i are all different, such that no two qubits are near a 
degeneracy, the interaction (|l5|) has only weak effects. In 
this case, the eigenstates of, say, a two-qubit system are 
approximately | J. J. ) , | J. f ) , ||4) and Iff), separated by en- 
ergies larger than El- Hence, the effect of the coupling 
is weak. If, however, a pair of these states is degener- 
ate, even a weak coupling lifts the degeneracy, changing 
the eigenstates drastically. For example, if V\ = V2, the 
states ||t) an d If I) are degenerate. In this case the cor- 
rect eigenstates are: ^(|IT) + |U» and ^(UT)-ITD) 
with energy splitting El between them. 

With the coupling ( |l5| ) we are able to perform two-bit 
operations and create, e.g., entangled states. In the idle 
state we bias all qubits at different voltages. Then we 
suddenly switch the voltages of two selected qubits to be 
equal, bringing those two qubits for a time At into reso- 
nance, and then we switch back. The result is a two-bit 
operation, which is a rotation in the subspace spanned 
by lit), ITI): 



U 2 -ut{At) 



cos (3 i sin (3 
i sin (3 cos (3 



(16) 
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where [3 = E L At/2h. The states |U),|TT) merely ac- 
quire phase factors, relative to what they would acquire 
in the idle state. For all other qubits the interaction re- 
mains a small perturbation. 

Similar to the situation with single-bit gates, the many- 
qubit states have time-dependent phase factors and the 
result of consecutive 2-bit gates depends on the waiting 
time between them. Again a transformation to the in- 
teraction representation makes this time dependence ex- 
plicit. One can show that any two-qubit operation can, in 
principle, be performed exactly. However, the complex- 
ity of such exact manipulation grows with the number of 
qubits. The improved design proposed in [fl9[ allows one 
to switch on and off the coupling between the qubits and 
to perform exact two-qubit operations in a simple way. 

The two-bit gates (|l6|) together with all one-bit gates, 
i.e. spin rotations (^|) and the phase shifts, constitute a 
universal set: They are sufficient for all manipulations 
required for quantum computation. 



tions of the design of the qubits and their coupling may 
be useful. We mention here one which has been described 
in Ref. jl9| . There we discuss a system where the Joseph- 
son coupling can be tuned to zero as well. This can 
be achieved by replacing each Josephson junction by a 
SQUID with two junctions, and controlling the effective 
coupling by an applied magnetic flux. While this design is 
more complicated to fabricate, it has considerable advan- 
tages: (i) In the idle state the Hamiltonian can be tuned 
to zero (E c h(V q b) = -Ej.cff = 0), which removes the prob- 
lem that the unitary transformations not only depend 
on the time span of their duration but also on the time 
to when they are performed, (ii) The two-bit coupling 
is turned on only for those two qubits which have both 
^j.eff 7^ 0. (iii) With only -E c h(V q b) or -Ej^ff non-zero, the 
unitary transformations depend in a simple way on the 
time integral of the corresponding Hamiltonian. Hence 
ramping the energies produces simple, well-defined gates. 
As a result of these improvements the manipulations can 
be performed with much higher accuracy. 



C. Extensions and discussion 



To check the experimental feasibility of the present 
proposal it is necessary to estimate the time-span At 
of the voltage pulses needed for typical single-bit opera- 
tions. We note that a reasonable value of Ej is of the or- 
der of Ej/Icb = 0.1-1K. It cannot be chosen much lower, 
since the condition k^T < Ej should be satisfied, and 
it should not be much larger since this would increase 
the technical difficulties associated with the time con- 
trol. The corresponding time-span is nevertheless short, 
At w h/Ej = 10-100ps, and is difficult to control in an 
experiment. 

There exist, however, alternatives. On one hand, a 
controlled ramping of the coupling energy provides a well 
defined, though non-trivial single-bit operation. Since in 
general, any 1-bit gate can be performed with proper 
choice of 3 controlled parameters, a universal set of gates 
can be produced in this way. Another possibility is to fol- 
low the established procedures of spin resonance experi- 
ments. E.g. a coherent spin rotation can be performed as 
follows: The system is moved adiabatically to the degen- 
eracy point. Then an ac voltage with frequency Ej/h is 
applied. Finally, the system is moved adiabatically back 
to the idle point. The time- width of the ac-pulse needed, 
e.g., for a total spin flip depends on the ac-amplitude, 
therefore it can be chosen much longer then h/Ej. Un- 
fortunately, in comparison to the sudden switching, the 
slow adiabatic approach allows only a reduced number of 
operations during the phase coherence time. 

Also the two-bit gates, instead of application of short 
voltage pulses, can be performed by moving the system 
adiabatically to a degeneracy point (say V\ = V2), and 
then applying an ac voltage pulse in the antisymmetric 
channel {V\ — V2) oc exp(i-Ef,t). 

To improve the performance of the device generaliza- 



III. CIRCUIT EFFECTS: DISSIPATION AND 
DEPHASING 



A. Johnson-Nyquist noise of the gate voltage circuit 



The idealized picture outlined above has to be ex- 
tended to account for the possible dissipation mecha- 
nisms causing decoherence and relaxation. We focus on 
the dissipation and fluctuations which originate from the 
circuit of the voltage sources. In Fig. || the equivalent cir- 
cuit of a qubit coupled to an impedance Z(w) is shown. 
The latter is characterized by intrinsic voltage fluctua- 
tions (between its terminals when disconnected from the 
circuit), with spectrum 



00 -I 

(6V8V) U = I dte Mt -(5V(t)SV(Q) + 6V(0)SV{t)) 
-00 ^ 



= Rc{Z(lu)}Tllu coth 



2k B T 



(17) 



When Z(u>) is embedded in a circuit, similar to that of 
Fig. H but with Ej = 0, the voltage fluctuations between 
the terminals of Z{ui) are characterized by a modified 
spectrum: 



(5VSV) U = Re{Z t (oj)}hcu coth 



TlLO 

2k B T 



(18) 



Here Zt(u>) = \iwCt + Z 1 is the total impedance 

between the terminals of Z(ui). 
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B. Relaxation and dephasing rates 



E, 



Z(co) 



FIG. 3. Qubit with electromagnetic environment. 

Following Caldeira and Leggett |||] we model the dis- 
sipative element Z(uj) by a bath of harmonic oscillators 
described by the Hamiltonian 



bath 



E 



i 



m a u>„ 



2m a 



(19) 



It is assumed that the voltage between the terminals of 
Z(u>) is given by SV = ^a%a and the spectral func- 

A 2 

tion J{uo) = § 2J Q " S(ui — io a ) is chosen to repro- 
duce the fluctuations spectrum (|l7|), J{uj) = uiRe{Z(u)}. 
Then, we embed the element Z(uj) in the qubit's circuit 
and we use the Kirchhof constraints to derive the Hamil- 
tonian of the whole system. Taking into account the 
possibility of a time-dependent external voltage V q b(t), 
we arrive at the resulting Hamiltonian 



H 



(2en) 2 



G 



2C*j 
2en 

~c7 



COS © + -ffbath 



(20) 



where the tilde-marking of the bath-related quantities 
reflects the fact that the bath has been modified and 
its spectral function corresponds now to the fluctuations 
spectrum (p"8|), J(lu) = LoRe{Z t {u))} . Shifting the ori- 
gins of all the oscillators x a — > x a — (A a /m Q d)„)Ctlqb,0) 
where V^b.o is the time independent part of the external 
voltage, we get a more familiar form of the Hamiltonian: 



ti — — -frj COS W + -Hbath 



c t 



2C 3 
2en 

c7 + 



SV, b (t) 



(21) 



To be specific we will concentrate in the following on the 
fluctuations due to an Ohmic resistor Z(oj) — R in the 
bias voltage circuit. 



In general, the environment has two effects: inelastic 
energy relaxation and dephasing, both characterized by 
their respective time scales, Fr and TV 1 . To illustrate 
this point we consider a qubit prepared in a superposi- 
tion of eigenstates a|0) + 6|1) with non-vanishing coef- 
ficients a and b. This corresponds to an initial density 

matrix a = (N' ft). In this case, one question is how 

fast the diagonal elements of the density matrix relax to 
their thermal equilibrium values. This relaxation is de- 
termined by T nl . The second question is how fast the 
off-diagonal elements vanish, which is governed by the 
rate T^. In general, the two rates are not equal. 

Both inelastic transitions and dephasing were ad- 
dressed in the context of spin-boson models, in particu- 
lar for a two-level system with purely Ohmic dissipation 
p9| , j30|] . Our system ( |2l| ) reduces to a spin-boson model 
in the limit where only two charge states need to be con- 
sidered, with a harmonic oscillator bath coupling to a z . 
It still differs from an Ohmic model due to the frequency 
dependence in the function Z t (uj). However, for the rel- 
evant frequencies, u) <C l/RCt, the physics is dominated 
by the resistor. Hence we can take over the results from 
the model calculations of Refs. |^9|,[30| for purely Ohmic 
environments. 

The strength of the fluctuations and, hence, the relax- 
ation and dephasing effects are controlled by the resis- 
tance of the circuit. It turns out that it has to be mea- 
sured in units of the quantum resistance Rk = h/e 2 ~ 
25.8kf2. Hence for low circuit resistances in the range 
of R « 5017 we can expect a weak effect. Further- 
more, in our system the effect of fluctuation is reduced 
due to the weak capacitive coupling of the qubit to the 
circuit. Indeed, as is apparent from ( pTh the coupling 
of ^2 a X a x a to the qubit's charge n involves the ratio 
C/(Cj + C) = Ct/Cj < 1. Thus, the relevant parame- 
ter characterizing the effect of the voltage fluctuations is 
(R/R K ) {CJCjf. 

In Ref. [£9| only the unbiased case of the spin-boson 
model, corresponding to the qubit at the degeneracy 
point n q b = 1/2, has been analyzed rigorously in the limit 
of low dissipation (the only limit relevant for the quan- 
tum computation). The biased case of the spin-boson 
model at low dissipation (i.e. the qubit out of the degen- 
eracy) was treated in Ref. jj(J. The two rates, expressed 
in terms of the mixing angle introduced in (|5|), are 



R ( C 

F in = sin 2 r) 47r ^r; ( ) -zr coth 



AE 

T 



AE 
2k B T ' 



and 



cos r\ 8ir 



R 
Rk 



C, 



(22) 



(23) 



We note that out of the degeneracy the dephasing rate 
acquires a component proportional to the temperature. 
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It should also be noted that the relaxation rate ( p2| ) may 
be obtained in a much simpler way JlG] using the so called 
"P(£)" theory discussed in Refs. |^f§4|. 

Thus, the diagonal and off-diagonal elements in the 
basis of the qubit's eigenstates relax with different rates 
described above, while the off-diagonal elements carry 
also an oscillating phase factor. The last is translated 
in the charge basis to the 'coherent charge oscillations', 
which are observed in the quantity (a z (t)) . In the absence 
of dissipation this quantity oscillates coherently around 
some average value. With dissipation, this average value 
relaxes to equilibrium with rate Tin, while the oscillations 
decay with rate |30j . 

The factors sin 2 r\ and cos 2 rj in Eqs. ( p^ , p3| ) indicate 
the nature of different contributions to the relaxation and 
dephasing rates. In the basis of the eigenstates (||) the 
qubit's part of the Hamiltonian ( pl| ) is diagonal. How- 
ever, the circuit electromagnetic fluctuations couple to 
the charge. Hence, in the notation introduced after (||) 
we have 



H 



9i 

Cj 



2e a z \ a i. 



= -^-2e (cos?7p z + s\wq p x )^^\ a x a . (24) 



Since the first term commutes with the qubit's Hamil- 
tonian, only the second one in J24| ) contributes to the 
transitions between the eigenstates. This explains the 
factor sin 2 r\ in (53) . The first term in ( jgjj ) , while ineffec- 
tive for transitions, makes the energy difference between 
the qubit's levels fluctuate. Therefore, the off-diagonal 
elements of the density matrix acquire a random phase. 
The last process is "pure" dephasing. It occurs even when 
transitions are suppressed (e.g. when Ej — 0). This par- 
ticular case (Ej = 0, 8Vqb(t) = 0) is quite simple and we 
can reproduce Eq. (|23| ) in an exact calculation, examin- 
ing also the zero temperature limit. To do this, we note 
that in the two state approximation the off-diagonal ele- 
ment of the reduced density matrix is given by <7o.i(f) = 
(exp(— iQ(t))}. To calculate it we perform a canoni- 
cal transformation x' a = x a — (C t /Cj)(A Q /m Q ,(D 2 ) 2en, 
6' = e + (27r/$ )(Ct/Cj) Y, a (^a/m a Ql)p a , after which 
the Hamiltonian (|2l|) becomes: 



H 



(2en-Cy qb , ) 2 
2(C + Cj) 

Ej COS 6 - — — > -r—^Vo, 



Et _/ 2en 
Kx' a + — 



(25) 



The Josephson and the last terms are zero in 
our case and the evolution of {0',ro} is decoupled 
from the bath. In the two state approximation 
(n 2 = n) the equations of motion read: dn/dt = 



0, d[exp(— iQ')]/dt = — iE c h(^jb,o) exp(— where 
E c h(..) was introduced in Eq. ([|). Thus, introducing 
$ = (2V$o)(C t /Cj)^ Q (A Q /m Q cD 2 )p Q , we get 

<7 0)1 (t) = <exp(i [$(t) - 0'(t)])) 

= (exp(i$(t)) exp(-i [9'(0) + ^ ch (y qb ,o) t])> 

= ( T ,i(D) e "' Jicl(y,b '° )t (e** (t) e- l$(0) ) . (26) 

The correlator {e 1 ®^ e -1 *^ -*) was studied extensively by 
many authors pl|-|34|]. It is equal to exp(4(Ct/Cj) 2 J(t)), 
where 



J(t) = 2 
x coth 



duj ReZ t (uj) 
o w Rk 



2k B T 



cos(o;t) — 1] — i sin(wt) 



(27) 



Since ReZ t (uj) = i?/(l + i? 2 C t 2 o; 2 ), we may roughly sub- 
stitute the bath's spectrum by a purely Ohmic one with 
a cut-off set at lo c = uirc — \/{RCt). Then, at non- 
zero temperature, ReJ(i) — (27rfceT '/Ti)(R/ R^jt for 
t > h/2kBT, and we reproduce Eq. ( ^3|) in the limit Ej — 
0. At zero temperature ReJ(i) s» —(2R/Rk) \og(uj c t) 
for t > l/u> c . Thus, the fact that the "pure" dephas- 
ing rate vanishes at zero temperature does not mean 
that there is no dephasing at all. The decay of the off- 
diagonal element of the reduced density matrix is just 
non-exponential but rather algebraic (see Refs. [p5 36|). 

We conclude that the dephasing of the initial quantum 
state of the qubit is caused by two different processes: 
the dissipative transitions between the eigenlevels, and 
the fluctuations of the energy difference between the lev- 
els. The off-diagonal elements of the density matrix are 
suppressed by both of them. While the first process sur- 
vives at low temperature, the second one produces the 
dephasing rate proportional to the temperature. 

These results should be taken into account when select- 
ing the optimum idle state for the qubit. If the temper- 
ature is low, k-gT -C Ej, the best choice is obviously far 
from the degeneracy point (sinry <C 1). Also at higher 
temperatures this regime has a longer coherence time 
than the degeneracy point, although only by a numer- 
ical factor. In all cases it is pure dephasing, rather than 
inelastic transitions which limits the coherence at the idle 
point. 



C. Discussion and extensions 

In summary, the dephasing rate is small if the effective 
resistance of the circuit is low compared to the quantum 
resistance, Rk- Furthermore, a low gate capacitance C 
reduces the coupling of the qubit to the environment. 
Hence, with suitable parameters (R < 5017, C/Cj < 0.1) 
at low temperatures the number of operations which can 
be performed before the environment destroys the phase 
coherence may be as large as 10 3 -10 4 . The value of the 
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phase coherence time in this case is of the order of 10- 
100 ns. Coherence times of this order of magnitude have 
been observed in experiments on quantum dots pffl. 

Longer phase coherence times would be achieved if the 
coupling of the qubit to the noisy environment could be 
further reduced, without weakening the coupling to the 
common inductor L which provides the 2-bit coupling. 
The model discussed above is not ideal since the Ohmic 
resistor with Johnson-Nyquist noise is assumed to be in 
the immediate vicinity of the qubit, reduced in its ef- 
fect only by a low gate capacitor. A suitable design with 
a combination of superconducting leads and filters can 
drastically improve the situation. 

If the physical dephasing effects are reduced, the more 
serious will be the errors related to imperfections in the 
time control or imprecise parameters and coupling ener- 
gies, e.g., non- vanishing 2-bit couplings during the idle 
periods. If the combination of all these effects is suffi- 
ciently reduced, allowing for 10 4 -10 5 coherent manipu- 
lations steps, then eventually the remaining errors can 
be corrected by suitable codes |38|| . It appears from our 
analysis that this goal can be reached with Josephson 
junction qubits. 



IV. MEASURING THE STATE OF THE QUBIT 



The model for a single-electron transistor 
attached to the qubit 



The read-out of the state of the qubit requires a quan- 
tum measurement process. Since the relevant quantum 
degree of freedom is the charge of the qubit island, the 
natural choice of measurement device is a single-electron 
transistor (SET) . This system is shown in Fig. |^. The left 
part is the qubit, with state characterized by the number 
of extra Cooper pairs, n, on the island, and controlled 
by its gate voltage, V^t,- The right part shows a normal 
island between two normal leads, which form the SET. 
It's charging state is characterized by the number of ex- 
tra single-electron charges on the middle island, N . It 
is controlled by gate and transport voltages, V g and Vtr, 
and further, due to the capacitive coupling to the qubit, 
by the state of the latter. A similar setup has been stud- 
ied in the experiments of Refs. |i"T| with the purpose to 
demonstrate that the ground state of a single Cooper pair 
box is a coherent superposition of different charge states. 
We will discuss the relation of these experiments to our 
proposal below. 

During the quantum manipulations of the qubit the 
transport voltage Vt r across the SET transistor is kept 
zero and the gate voltage of the SET is chosen to tune 
the island away from degeneracy points. Therefore no 
dissipative currents flow in the system, and the transis- 
tor merely modifies the capacitances of the system. To 
perform a measurement one tunes the SET by V g to the 



vicinity of its degeneracy point and applies a small trans- 
port voltage V tI . The resulting normal current through 
the transistor depends on the charge configuration of the 
qubit, since different charge states induce different volt- 
ages on the middle island of the SET transistor. In order 
to investigate whether the dissipative current through 
the SET transistor allows us to resolve different quan- 
tum states of the qubit, we have to discuss various noise 
factors, including the shot noise associated with the tun- 
neling current and the measurement induced transitions 
between the states of the qubit. For this purpose we 
analyze the time evolution of the density matrix of the 
combined system. We find that for suitable parameters, 
which can be realized experimentally, the dephasing by 
the passive SET is weak. When the transport voltage is 
turned on the dephasing is fast, and the current through 
the transistor — after a transient period — provides a 
measure of the state of the qubit. At still longer times 
the dynamics of the SET destroys the information of the 
quantum state to be measured. 



Vtr 
2 



Ej 
Cj 



m / C„ 



CintiH 

n "~ I I N 



C 



Vtr 
2 



FIG. 4. The circuit consisting of a qubit plus a SET tran- 
sistor used as a measuring device. 

The Hamiltonian of the combined system 



H — i? c h + Hi, + Hn + Hi 
+ Hi + H T 



(28) 



contains the charging energy, the terms describing the 
microscopic degrees of freedom of the metal islands 
and electrodes, and the tunneling terms, including the 
Josephson coupling. The charging term is a quadratic 
form in the variables n and N. I.e. for the system shown 
in Fig. [| it is 

H ch (n, N, V n , V N ) = 4E qh n 2 + E SET N 2 + 2E int nN 

+2enV n + eNV N . (29) 
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The charging energy scales -E q b, -Eset and E{ ai are de- 
termined by the capacitances between all the islands. In- 
troducing 

A = {C + C, 1 )(C g + C hlt + 2C / ) + C int (C g + 2C') « 2C 3 C , 

we can write them as 

E qh = e 2 (C g + CW + 2C , )/2A » e 2 /2Cj , 
^set = e 2 (C + C int + Cj)/2A « e 2 /(4C') , 
E int = e 2 C int /A « e 2 C int /(2CjC) . (30) 

Here we have assumed that the two junctions of the SET 
have equal capacitances C , and the approximate results 
refer to the limit C, C- mi ,C g «C'« Cj , which we con- 
sider useful (see below). The effective gate voltages V n 
and Vn depend in general on all voltages V q b, V e , and 
the voltages applied to both electrodes of the SET. How- 
ever, for a symmetric setup (equal junction capacitances) 
and symmetric distribution of the transport voltage, Vt T , 
between both electrodes of the SET (as shown in Fig. 
V n and Vjv are controlled only by the two gate voltages, 



V N = V 



C g (C + C\nt + Cj) . t r Ci n tC 



A 



+ V qh - 



A 



r„ = v c&* + Vqh (c g + c in ^+2C')c , Tt 



A 



The microscopic terms describe noninteracting elec- 
trons in the two leads and on the middle island of the 
SET transistor 



E«t4 CT (r = L,R,I) 

A: a - 



(32) 



The index a labels transverse channels including the spin, 
while k labels the wave vector within one channel. 

Similar terms describe the electrode and island of the 
qubit; however, for the superconducting non-dissipative 
element the microscopic degrees of freedom can be in- 
tegrated out 39 1, resulting in the "macroscopic" quan- 
tum description presented in Sect. II. In this limit the 
tunneling terms reduce to the Josephson coupling Hj = 
— -EjcosO, expressed in a collective variable describ- 
ing the coherent transfer of Cooper pairs in the qubit, 
e i@ \n) = |n+l). 

The normal-electron tunneling in the SET transistor is 
described by the standard tunneling Hamiltonian, which 
couples the microscopic degrees of freedom, 



L-i 1 kk' a v ka L 'k' a c 



-i(f> 



kk' 



k'k" 



(33) 



To make the charge transfer explicit, ( 

a ±id> QT1 J „±i 



u displays two 

macroscopic" operators, e =ci<? ' and e ±l ^T The first one 
describes changes of the charge on the transistor island 



due to the tunneling: e l<l} \N) = \N + 1). It may be 



treated as an independent degree of freedom if the total 
number of electrons on the island is large. We further 
include the operator e ±l ^ which describes the changes of 
the charge in the right lead. It acts on m, the number 
of electrons which have tunneled through the SET tran- 
sistor, e l ^\m) = \m + 1). Since the chemical potential 
of the right lead is controlled, m does not appear in any 
charging part of the Hamiltonian. However, we have to 
keep track of it, since it is the measured quantity, related 
to the current through the SET transistor. 

In equilibrium, i.e. when the SET is kept in the state 
N = 0, the qubit 's dynamics is described by the same 
Hamiltonian as discussed in previous sections, i? q b = 
4£; qb (?i - n qh ) 2 - EjcosQ, where n qb = -eV n /4E qh . 
We recall that in the limit where only the lowest en- 
ergy charge states n = and n = 1 are relevant, it re- 
duces to a two state quantum system. In the basis of 
eigenstates (||), |0) and |1), which are expressed in terms 
of the mixing angle ?y, where tan7y = Ej/ E cn (n q h), it 
becomes H q \> — ^AE V p z , where AE n = Ej/smr] and 



qbj 



4£ q b(l — 2n q b). In this basis the number 



E ch (n, 

operator n becomes non-diagonal, 

n = ^(1 + <t z ) = -(1 + cosi]p z + sinr]p x ) . (34) 

For the following discussion we choose n q b away from 
the degeneracy point, which combined with Ej <§C -E q b 
implies tan 77 <C 1. 



The mixed term in (29) provides the interaction Hamil- 
tonian. With n given by (|i]) it becomes 

Hint = E iDt Na z = E int N(cosr]p z + siwq p x ) , (35) 

plus an extra term, Ei nt N, which together with further 
terms is collected in the Hamiltonian of the SET transis- 
tor, 



-ffsET = E$et(N — Nset) 2 + + + H\ + Ht 



(36) 



The transistor's gate charge became ./Vset = — (eVjv 
E- m t)/2EsET- The total Hamiltonian thus reads 



H = H, 



qb 



SET 



H, 



(37) 



The total system composed of qubit and SET is de- 
scribed by a total density matrix p(t), which we can re- 
duce, by taking a trace over the microscopic states of 
the left and right leads and of the island, to a(t) = 
TrL,R,i{/3(£)}. In general, this reduced density matrix 
a(i,i'; N, N';m,m') is a matrix in the index i, which 
stand for the quantum states of the qubit |0) or |1), in N, 
and in m. In the following we will assume that initially 
— as a result of previous quantum manipulations — the 
qubit is prepared in the quantum state a|0) + and 
at time t = we switch on a transport voltage to the 
single-electron transistor. We then proceed to further re- 
duce the density matrix in two ways which provide com- 
plementary information about the measuring process. 
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The first, widely used procedure [£3| is to trace over 
N and to. This yields a reduced density matrix of the 
qubit Oi.j = J2n m^"(*' 3> N \ N;m,m). Just before the 
measurement, it is in the state 



*(0) 



a*b 



ab* 
\b\ 2 



(38) 



The questions then arise how fast the off-diagonal ele- 
ments of o"jj vanish, i.e. how fast is the dephasing, and 
how fast the diagonal elements change their values af- 
ter the SET is switched to the dissipative state. These 
questions are the sa me a s in our discussion of fluctua- 
tion effects in Section III. However, the coupling to the 
dissipative SET in general shortens these times. This 
description is enough when one is interested in the quan- 
tum properties of the measured system, i.e. the qubit 
only, and the measuring device serves merely as a source 
of dephasing [fjj"|-p4f . It does not tell us much about the 
quantity measured in an experiment, namely the current 
flowing through the SET. 

The second procedure, pursued in the following, is to 
evaluate the probability distribution of the number of 
electrons to which tunnel through the SET during time 
t 



P(m,t) 



N,i 



a(i, i; N, N; m, m)(t) 



(39) 



This distribution provides the information about the 
experimentally accessible quantity during the measure- 
ment process. At t = no electrons have tunneled, so 
P(m, 0) = 5 m fi. Then the peak of the distribution moves 
in positive m-direction and, simultaneously, it widens due 
to shot noise. Since two states of the qubit correspond to 
different tunneling currents, and hence shift velocities in 
m-direction, one may hope that after some time the peak 
splits into two. If after sufficient separation of the two 
peaks their weights are still close to \a\ 2 and \b\ 2 , a good 
quantum measurement can be performed by measuring 
to. After a longer time further processes destroy this 
idealized picture. The two peaks transform into a broad 
plateau, since transitions between the qubit's states are 
induced by the measurement. Therefore, one should find 
an optimum time for the measurement, such that, on one 
hand, the two peaks are separate and, on the other hand, 
the induced transitions have not yet influenced the pro- 
cess. In order to describe this we have to analyze the time 
evolution of the reduced density matrix quantitatively. 



B. Quantitative description of the measurement 
process 

The time evolution of the density matrix leads to 
Bloch-type master equations with coherent terms. Ex- 
amples have recently been analyzed in contexts similar 
to the present [|ol|4^J2|]. In Rcf. Q a diagrammatic 



technique has been developed which provides a formally 
exact master equation as an expansion in the tunneling 
term Ht, while all other terms constitute the zeroth or- 
der Hamiltonian Hq = H — , which is treated exactly. 
The time evolution of the reduced density matrix is given 
by &(t) — <r(0)n(0, t). The propagator H(t' , t) can be ex- 
pressed in a diagrammatic form and finally summed up 
in a way reminiscent of a Dyson equation. Examples are 
shown in Fig. [5| In contrasts to ordinary many-body ex- 
pansions, since the time dependence of the density matrix 
is described by a forward and a backward time-evolution 
operator, there are two propagators, which are repre- 
sented by two horizontal lines (Keldysh contour). The 
two bare lines describe the coherent time evolution of 
the system. They are coupled due to the tunneling in 
the SET, which is treated as a perturbation. The sum of 
all distinct transitions defines a 'self-energy' diagram £. 
Below we will present the rules how to calculate £ and 
present a suitable approximate form. The Dyson equa- 
tion is equivalent to a (generalized) master equation for 
the density matrix, which reads 



da(t) 
dt 



&(t),H ] = / dt'a(t')T,{t-t') 



(40) 



In the present problem the density matrix is a matrix 



~i ,N ,m 



i, N, m 



in all three indices 



a(i,i';N,N';m,m') = a 
N, and to, and the (generalized) transition rates due 
to single-electron tunneling processes (in general of ar- 
bitrary order), S 1 ' : ™ ~~ ,' m (t — t 1 ), connect these 
diagonal and off-diagonal states. 



n 



n 



(0) 



n 



n 



(0) 



FIG. 5. The Dyson-type equation governing the time evo- 
lution of the density matrix. It is equivalent to the generalized 
master equation mw. The 'self energy' diagrams E describes 
the transitions due to tunneling in the SET transistor. 

The transition rates can be calculated diagrammati- 
cally in the framework of the real-time Keldysh contour 
technique. We briefly review the rules for their evalua- 
tion; for more details including the discussion of higher 
order diagrams we refer to Ref. pOj . Typical diagrams, 
which will be analyzed below, are displayed in Fig. ^ 
and[7j A gain the horizontal lines describe the time evo- 
lution of the system governed by the zeroth order Hamil- 
tonian Hq- Their properties will be discussed below. The 
directed dashed lines stand for tunneling processes, in the 
example considered the tunneling takes place in the left 
junction. According to the rules the dashed lines con- 
tribute the following factor to the self-energy S, 



irk B T\ 2 exp[±i^ L (*-*') 



sinh 2 (t-t'± iS)] 



(41) 



11 



where = h / (Air 2 e 2 Rt ,l) is the dimensionless tunnel- 
ing conductance, /il is the electro-chemical potential of 
the left lead, and <5 _1 is the high-frequency cut-off, which 
is at most of order of the Fermi energy. The sign of 
the infinitesimal term iS depends on the time-direction 
of the dashed line. It is negative if the direction of the 
line with respect to the Keldysh contour coincides with 
its direction with respect to the absolute time (from left 
to right), and positive otherwise. For example the right 
part of Fig. || should carry a minus sign, while the left 
part carries a plus sign. Furthermore, the sign in front of 
ifjb^lt — t') is negative (positive), if the line goes forward 
(backward) with respect to the absolute time. The first 
order diagrams are multiplied by (—1) if the dashed line 
connects two points on different branches of the Keldysh 
contour. 



N,m,j'|- 

N,m,j'-»N+l,mj 
N,m,i'— N+l,m,i 

N,m,i':'- 



N+l,mj 


N,m,j' 








t + 


t' 




\ L 


N+l,m,i 


N,m,i" 









N+l,mj 



FIG. 6. Example of a 'self energy' diagram for an "in" 
rate. 



N.mj'j- 



N,m,j'-^N,m,j 



N,m,i'— N,m,i 



t' 



N,m,j N,mj' 




v *7 


t + t' 

N,m,i N,m,i' 


L 







N,m,j 
t 

N,m,i 



FIG. 7. Example of a 'self energy' diagram for an "out" 
rate. 

The horizontal lines describe the time-evolution of the 
system between tunneling processes. For an isolated 
transistor island they reduce to simple exponential fac- 
tors e ±: k E ( t ~ t \ depending on the charging energy of 
the system. In the present case, however, where the 
island is coupled to the qubit we have to account for 
the nontrivial time evolution of the latter. For instance, 
the upper line in the left part of Fig. [| corresponds to 
(iV, j\e~ : k H °^ t ~ t '\N,j'}, while the lower line corresponds 
to (N + l,i'\ef: Ho( - t ~ t ")\N + l,i). 

In principle the density matrix is an arbitrary non- 
diagonal matrix in all three indices i, N, and m. But, 
as has been shown in Ref. |40[| , a closed set of equa- 
tions can be derived, describing the time evolution of 
the system, which involves only the diagonal elements 
in N. The same is true for the matrix structure in m. 
Therefore, we need to consider only the following ele- 
ments of the density matrix cr^'m- Accordingly, of all 
the transition rates we need to calculate only those be- 
tween the corresponding elements of the density matrix, 



i.e. £ , 



i' ,N' 1 Tn t — >i,N ,m 



j' ,N' ,m' — >j.N,m 



(At). In the present problem we fur- 
ther can assume that the tunneling conductance of the 
SET is low compared to the inverse quantum resistance. 



In this case, lowest order perturbation theory in the single 
electron tunneling, describing 'sequential tunneling pro- 
cesses', is sufficient. The diagrams for £ can be split into 
two classes, depending on whether they provide expres- 
sions for off-diagonal (N' ^ N) or diagonal (N' — N) in 
TV elements of £. In analogy to the scattering integrals 
in the Boltzmann equation these can be labeled "in" and 
"out" terms, in the sense that they describe the increase 
or decrease of a given element Sj'jy ™ of the density ma- 
trix due to transitions from or to other TV-states. Exam- 
ples for the "in" and "out" terms are shown in Fig. |^ and 
Fig. 0, respectively. 

We now are ready to evaluate the rates in Fig. || and 
Fig. 0. As an example we consider an "in" tunneling 
process in the left junction: 



-if ,N,m—>j,N+l,m 
•V ,N,m—*i,N+l,m 



(At) 



exp 



i ( fpN 



W» +1 )At 



exp 



sinh 2 [=^(At 



■iS)] 



(42) 



pN-\ 



w 



jV+1' 
N 



At 



j j 



sinh 



2 [irk B T 



(At - iS)] 



where E 



(E SET (N 1 ) -// L JVi) - (E SET (N 2 ) -fx L N 2 ) 



is the usual Coulomb energy gain for the tunneling in the 
left junction in the absence of the qubit, and 



W$ = fl&(JVi) ® 1 - 1 ® H qh (N 2 ) 



(43) 



provide corrections to the energy gain sensitive to the 
state of the qubit. Here H q ^(N) is the iV-th block of 
the Hamiltonian i/ q b + -ffmt (note that 7tf q b and H lnt are 
block-diagonal with respect to N). The indices j',j and 
i', i relate to the left and right side of the tensor product 
in (^3|) correspondingly. 

The form of the master equation (^0|) suggests the use 
of the Laplace transform, after which the last term in 
( fto| ) becomes Y,(s)a(s). We Laplace transform (|42] ) in the 
regime hs, \ W$ +1 \, \ W^ +1 \ <C E% +1 , i.e. we assume the 
density matrix a to change slowly on the time scale given 
by fr/E$ +1 . This assumption should be verified later for 
self-consistency. These inequalities also mean that we 
choose the operation regime of the SET far enough from 
the Coulomb threshold. Therefore, the tunneling is either 
"allowed" for both states of the qubit or it is "blocked" 
for both of them. At low temperatures (fceT <§C E$ +1 ) 

and for E$ +1 S -C Ti we obtain: 



^N,m,j f —*-N-\-l,m } j/ \ 
^jV.m.i'^AT+l,™,,* \ S ) 



-a L 0« +1 ) 



2E 



N 

N+l 



-a L D(E») 



(44) 
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where D{E™ +1 ) w 1 + 7 + ln(£$ +1 <5//i) and 7 w 0.58 is 
Euler's constant. The first term of (Q) is the standard 
Golden rule tunneling rate corresponding to the so-called 
orthodox theory of single-electron tunneling ^2|. The 
rate depends strongly on the charging energy difference, 
-Etv+h before and after the process, which in the present 
problem is modified according to the quantum state of 
the qubit (the W terms). At finite temperatures the step- 
function is replaced by ®(E) — > [1 — exp (— _E/fc B T)] _1 . 
We denote the full matrix of such rates T and we exten- 
sively exploit this matrix below. If the temperature is 
low and the applied transport voltage not too high, the 
leading tunneling process in the SET is the sequential 
tunneling involving only two adjacent charge states, say 
N = and N — 1. We concentrate here on this case (to 
avoid confusion with the states of the qubit we will keep 
using the notation N and N + 1). 

The second, logarithmically diverging part of ( (fif ) pro- 
duces a "coherent-like" term in the RHS of the master 
equation (|40|). These terms turn out to be unimportant 
in the first order of the perturbation theory. Indeed, for 
the left junction we obtain the following contribution to 
the RHS of <M): 



da i 



(45) 



the simplest form. The matrix structure in the m 
variable simplifies considerably if we perform a Fourier 
transform with respect to to, crfj(fc) = 
To shorten formulas we introduce A N 



B 



JY 



N 



N,n 
0,1 



and D 



N 



Akra 



Im a. 



N,n 
0,1 



Then, introducing a vector X(k) 



(A n ,A n +\B n ,B n +\C n ,C n +\D n ,D n + 1 ), we can 
rewrite (E^) as 



X(k) = -M(k)X(k) , 
n 



where M(k) is given by 



/ — Tlo e lk ^R0 









-r L1 e lk r R1 
r L1 -r R1 



Ej 


-E 3 


















-\Ej ^Ej 

V -ijsj 






\E, 



E C ) 
E 



-Tr 


N+l 





r L 



and 



(48) 



-F L e ik T R -E^ 





E 3 


-Ej 



-E 



e"T K 

~Tr J 



(49) 



where H qh = ±(H qh (N) + H qh (N + 1)) and D L is a 
matrix in N and to spaces. The eigenvalues of the 
matrix Dl are at most of order D(E$ +1 ). Neglecting 
terms of order ai,D(E^ +1 )E in t in Eq. (|45|) , we can re- 
place .Sqb by H - Our analysis shows that these neglected 
"coherent-like" terms do not change the results as long 
as q;l I ln(Ejy +1 8/h)\ <C 1. Similar analysis can be made 
for the right tunnel junction of the SET. 

Now we can transfer all the "coherent-like" terms into 
the LHS of the master equation, 



1 - a L D L - a R D R 



jeZ<7(/i 



df 



[&{t),H ] 



1 



(46) 



and multiply Eq. (|46|) from the left by (1 — «l-Dl — 
c^rDr)^ 1 ~ (1 + + «rZ)r) so that the corrections 

move back to the RHS. Since F is itself linear in aL, a R 
the corrections belong to the second order in a's (more 
accurately, they are small if a\ln(E^ +1 S/h)\ <C 1 for 
both junctions). Thus, we drop the "coherent" correc- 
tions and arrive at the final form of the master equation 
which we use below: 



(47) 



To rewrite Eq. @ in a matrix form we have to 
choose a particular basis for the qubit. We choose 
the charge basis since in this basis the matrix T has 



T L = 2na h [ ml - ^set(1 - 2N SET )} , 

r R = 2™ R [- MR + £ S et(1 - 2N SET )} , (50) 

r L0 = T L + AT L , 
r L1 = T L - AT L , 

r R0 ee r R + at r , 

r R1 ee r R - at r . (51) 

The shifts Ar are proportional to the interaction energy 
between the qubit and SET, 



Ar L ee 2na L E int , 
Ar R ee -2ira R E int . 



(52) 



As we will see, these shifts are responsible for the separa- 
tion of the peaks of P(m, t). Finally, the qubit 's charging 
energies are given by 



E^ +1 = E ch (n qh ) + 2E h 



(53) 



C. Analysis of the master equation in the zeroth 
order in Ej. Measurement and dephasing times. 

We are interested in the limit Ej,Ei nt <C B c h(^qb) 
In this case the system ([l8]) may be analyzed perturba- 
tively in _Ej . In zeroth order the system of equations ( plSj 
factorizes into three independent groups. The first one, 
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TiA N+1 



-T L0 A N + T R0 e* k A N+1 , 
Tr, n A N - IW* +1 , 



(54) 



has exponential solutions cx e luJt and eigenvalues 



ROj 



1± 



4r (e 



ik 



i) 



(55) 



Here the rate Tq/H of single-electron tunneling through 
the SET if the qubit is in the state |0) is given by 



To 



TloTro 



r 



LO 



RO 



(56) 



When k is small the first eigenvalue, Tuji ~ j(IYo + Tro), 
has a large imaginary part and the corresponding eigen- 
mode quickly dies out. The second eigenvalue hu>2 ~ 
Tok + iTyCok 2 is small, and the corresponding eigenmode, 
with A N+1 /A N = r L0 /r R0 , survives. This solution de- 
scribes a wave packet propagating with the group velocity 
Tq, which widens due to shot noise of the single electron 
tunneling, its width being given by ^J^-T^t. The factor 
c = (T 2 L0 + r R0 )/(r L0 + r R0 ) 2 (see e.g. Ref. @) varies 
between 1/2 in the symmetric situation (Flo = Tro) and 
1 in the extremely asymmetric case (Flo much larger or 
much smaller than Tro)- 

Analogously the second group of equations for B N and 
B N+1 describe a wave packet which moves in m-direction 
with the group velocity 



r T .ir 



Ll-L Rl 



Tli + Tri 



(57) 



and the width growing as y/fTit The two peaks corre- 
spond to the qubit in the states |0) and |1), respectively. 
They separate when their distance is larger than their 
widths, i.e. jr\T -Ti\t > ^/fT^i+ ^/fT\i. This means 
that after the time 



( y/ToCQ + yT 



iCl 



V |ro-r x | 



(58) 



which we denote as the measurement time, the process 
can constitute a quantum measurement. 

To learn about the dephasing we analyze the last four 
equations of (|4^) at k = 0, which is equivalent to a trace 
over to. These four equations may be recombined into 
two complex ones: 



P. „• rpN ±.N F » jV , F ±.N+1 

fi « -JV+1 _ • - JV+1 , r -W r ~N- 

dt 0,1 ~ ch 0,1 4 ~ R o- 1 



N+l 



(59) 



The analysis of this set shows that if \E^ +1 — E^[\ = 
2E[ nt <C (Fl + Tr) the imaginary parts of the eigen- 
values are luihu>i ~ (T^ + Tr) and Im?kj2 ~ 



4^i 2 nt r L rR/(r L +T R ) 3 . In the opposite limit 2£ int > 
(Fl + Tr) the imaginary parts are Imhuji w Tl and 
Im T1UJ2 ~ Tr. The first limit is physically more relevant 
(we assume parameters in this regime) , although the sec- 
ond one is also possible if the tunneling is very weak or 
the coupling between the qubit and the SET transistor 
is strong. In both limits the dephasing time, which is 
defined as the the longer of the two times, 



T <t> 



max{[Imii)i] , [Ima>2] 1 } 



(60) 



is parametrically different from the measurement time 
@. In the first limit, 2E int < (T L + T R ), it is 



T ( j ) = h 



(r L 



r R ) 3 



4r L r R i?? t 



(61) 



One can check that in the whole range of validity of our 
approach the measurement time exceeds the dephasing 
time, t ms > T0. This is consistent with the fact that 
a "good" quantum measurement should completely de- 
phase a quantum state. 

The reason for the difference between the measurement 
time and the dephasing time is the entanglement of the 
qubit's state with additional microscopic states of the 
SET, — states which cannot be characterized by the 
number of electrons which have tunneled through the 
transistor, to, only. Imagine that we know (with high 
probability) that during some time from the beginning 
of the measurement process exactly one electron has tun- 
neled through the SET whatever the qubit's state was. 
Does it mean that with that high probability there was 
no dephasing? The answer is no. The transport of elec- 
trons occurs via a real state of the island, N+l. The 
system may spend different times in this intermediate 
state, i.e. different phase shifts are acquired between the 
two states of the qubit. Since the time spent in the state 
N + 1 is actually random, some dephasing has occurred. 
Next we can ask where the information about this phase 
uncertainty is stored, i.e. which states has the qubit be- 
come entangled with. These may not be the states with 
different to or N since m is equal to one and N is again 
equal to zero irrespective of the state of the qubit. The 
only possibility are the microscopic states of the mid- 
dle island and/or leads, which were subject to different 
time evolutions during different histories of the tunneling 
process. To put it in the language of Ref. p4j], the ini- 
tial state of the system (a|0) + 6|1)) |x) \fn = 0) evolves 
into a|0) |xo) |™o) + Ixi) l TO i)i where |x) stands for 
the quantum state of the uncontrolled environment. One 
may imagine a situation when too = mi, but |xo) and 
|xi) are orthogonal. In this situation the dephasing has 
occurred but no measurement has been performed. 

Dephasing was also analyzed in Refs. |2l| Q. In these 
works a quantum point contact (QPC) measuring device 
was used only as a source of dephasing, i.e. the informa- 
tion on the current flowing in the measuring device was 
disregarded. Thus a distinction between measurement 
and dephasing was not made. However the expressions 
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for the dephasing time were given b y e xpressions simi- 
lar to (58). As it became clear later |4j|, the QPC does 
not involve an additional uncontrolled environment dur- 
ing the measurement process and, therefore, t ms = t^,. 
Indeed, the tunneling in the QPC does not occur via a 
real intermediate state and the discussion above does not 
apply. In this sense the QPC may be regarded as a 100% 
efficient measuring device. The additional environment 
in the SET, which, as mentioned, reduces the efficiency 
of the measuring device, plays, actually, a positive role 
in the quantum measurement, provided it dephases the 
state of the qubit only when the system is driven out of 
equilibrium. This is because the quicker dephasing sup- 
presses the transitions between the states of the qubit 
(such a suppression of the transitions due to a contin- 
uous observation of the quantum state of the system is 
called the Zeno effect), while the initial probabilities are 
preserved. It is worth noting that a SET might also be 
used as a 100% efficient device, if it was biased in the 
co-tunneling regime. However, this possibility has other 
drawbacks and hence is not considered here further. 



D. Higher orders in Ej, the mixing time 

Finally, we analyze what happens if we take into ac- 
count the mixing terms proportional to Ej in the system 
(fl8|). We consider k = and investigate the eigenvalues 
of the 8x8 matrix ( flit ) • Note that in the discussion above 
we have calculated all the eight eigenvalues for Ej = 
(the two eigenvalues of the complex system (^9|) are dou- 
bled when one considers it as a system of four real equa- 
tions). In the diagonal part there were two zeros, which 
corresponded to two conserved quantities (for k = 0), 
A N (0) + A N+1 (0) = <7 , and B N (0) + B N+1 (0) = <7 M . 
Six other eigenvalues were large compared to the ampli- 
tudes of the mixing terms. It is clear, that including 
the mixing, Ej ^ 0, changes only slightly the values 
of the six large eigenvalues. A more pronounced effect 
may be expected in the subspace of the two degener- 
ate eigenvectors with zero eigenvalues. It turns out that 
this degeneracy is lifted in second order of perturbation 
theory. One of the resulting eigenvalues remains zero. 
This corresponds to the conservation of the total trace 
A N {0) + A N+1 (0) + B N (0) + B N+1 {0) = 1. The second 
eigenvalue acquires now a small imaginary part and this 
gives the time scale of the mixing between the two states 
of the qubit. In the limit of (our) interest we obtain 



t: 



4El t E]r 



where T 



?iA£ 2 (A£ 2 + (r L + r R ) 2 ) 

=s rs, and AE w AE(rj) 



(62) 



E ch 



E. 



N+l 
ch 



Comparing Eq. (|62| ) and Eq. (|58j ) we see that both cases 
t ms < imix and t ms > i mix are possible. 

Let us analyze a concrete physical situation. We as- 
sume oil = or = a. We choose N$et far enough 
from the degeneracy point, which is N$et — 1/2, so 



that Tl < Tr and the related Coulomb blockade energy, 
E C b = ^set(1 — 2Aset), is of the order of E sct . To sat- 
isfy the conditions for the Golden Rule (see ( |44| ) and the 
discussion thereafter) we assume Eqb to be the largest 
energy scale of the system, Eqb 3> AE, and the chemi- 
cal potential of the left lead, /iL = Vtr/2, to exceed the 
Coulomb blockade energy by an amount of the order of 
Eqb- The transport voltage should not, however, exceed 
the value beyond which further charge states of the SET 
transistor, e.g. N + 2 and N — 1, become involved. Thus 
Vtr/2 < E sct (l + 27Vset) and A^et should be chosen far 
enough from zero as well. Thus we obtain 



t: 



2TraE SET E? nt E 2 j 



h(AE) 2 max[(AE) 2 , (27raS S ET) 2 



(63) 



The measurement time in the same regime is given ap- 
proximately by 



t: 



H EsET 



(64) 



For comparison we also give a rough value for the de- 
phasing time, which is short when the measurement is 
performed. We assume the regime in which Eq. @ is 
valid. Then 



1 



T 4> 



E int 



2nafi -©set 



(65) 



Note that the limit a — > is not allowed in (|6^), since 
the validity of ( |6l| ) eventually breaks down in this limit. 
Thus, 



oc 



EjE s 



ET 



AE max[Ai;, 2naE SET } 



(66) 



One recognizes two competing ratios here: Ej/AE, 
which is small, and i5sET/max[A£', 27rai?sET] , which is 
large. The condition t ms /t m i x <C 1, thus, imposes an ad- 
ditional constraint on the parameters of the system. On 
the other hand, for low conductance barriers in the SET 
(a < 1 but not too small), the dephasing time is always 
shorter than the measurement time. 



E. Discussion and Extensions 

Let us summarize what has been done thus far: To 
measure the quantum state of the qubit we attach a SET 
to the qubit and consider the time evolution of the whole 
system. We provide two complimentary descriptions of 
the measurement process. In the first we trace over the 
SET's degrees of freedom and obtain a reduced density 
matrix of the qubit. From the time evolution of the last 
we deduce the dephasing time and the mixing time, i. e. 
the relaxation times for the off-diagonal and the diago- 
nal elements of the reduced density matrix respectively. 
This description does not, however, contain information 
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about the measurable quantity, the current in the SET. 
Therefore, in a second description, we calculate the prob- 
ability distribution, P(m,t), of the number of electrons 
m which have passed through the SET during time t. 
We will evaluate P(m, t) numerically for parameters to 
be specified in the next section. The results, shown in 
Fig. ^| and ^, display the time evolution on the time 
scale of t ms , relevant for the measurement process. 

We observe that under appropriate conditions P(m, t) 
splits into two separate peaks, whose weights are given 
by the initial probabilities of the eigenstates of the qubit. 
The splitting time is the minimum time after which one 
can distinguish between the states of the qubit and we 
call it, therefore, the measurement time, t ms . The mea- 
surement time turns out to be longer than the dephasing 
time. This fact indicates that the state of the qubit gets 
entangled not only with the transport degree of freedom 
in the SET but also with other, microscopic degrees of 
freedom which we do not observe fisj . 

Note, that the knowledge of P(m, t) does not provide 
immediately the value of the current flowing in the SET 
at all times. Recently, attempts have been made to de- 
scribe the measurement process in terms of the current 
in the measuring device |45|-[47j]. However, further inves- 
tigation of this problem is definitely needed. 



V. DISCUSSION 



A. Choice of parameters 



which can be performed. With circuit parameters dis- 
cussed below and for not too many qubits these restric- 
tions are weaker than those provided by the dephasing. 
Moreover, we note that these restrictions are removed 
if the interaction is effectively turned off. This can be 
achieved on the "software" level, by the use of refocusing 
voltage pulses. On the "hardware" level this is achieved 
with the improved design suggested in Ref. JlSfl . This 
allows for larger numbers of qubits and longer computa- 
tions. 

As an example we suggest a system with the following 
parameters: 

(i) We choose junctions with the capacitance Cj = 
4- 10~ 16 F, corresponding to the charging energy (in tem- 
perature units) -Eqb ~ 2K, and a smaller gate capacitance 
C = 2.5 ■ 10~ 18 F to reduce the coupling to the environ- 
ment. Thus at the working temperature of T — 50mK 
the initial thermalization is assured. The superconduct- 
ing gap has to be slightly larger A > E qD . We further 
choose Ej = lOOmK, i.e. the time scale of one-qubit op- 
erations is Top = h/Ej ~ 7 ■ 10~ n s. 

(ii) We assume that the resistor in the gate voltages cir- 
cuit has R ~ 50S7. Voltage fluctuations limit the de- 
phasing time (^3|). We thus arrive at an estimate of the 

decoherence rate which allows for (r^T^) -1 ~ 8 • 10 5 
coherent operations for a single bit. 

(iii) To assure sufficiently fast 2-bit operations we choose 
L ~ Then, the decoherence time is (r^ri^) -1 ~ 
650 times longer than a 2-bit operation. In Table 
[| we present alternative sets of parameters, assuming 
throughout that the resistance is fixed as R = 50fL 



To demonstrate that the constraints on the circuit pa- 
rameters can be met by available technology, we summa- 
rize the constraints and suggest a suitable set. The neces- 
sary conditions arc: A > E qD ^> Ej^El^-qT. The tem- 
perature has to be chosen low to assure the initial ther- 
malization, k^T <C E qD and k^T <C hujLc, and to reduce 
the dephasing effects. A good choice is k^T ~ Ej/2 since 
further cooling would not decrease the dephasing rate at 
the degeneracy point much further. Then, the dephasing 
rates during operations (at the degeneracy point) and in 
the idle state (off degeneracy) are close to each other. 

Depending on the parameters chosen, the number of 
qubits N and the total computation time are restricted. 
The dephasing times limits the number of operations to 
Top -C rr 1 , T^ 1 . This is the only restriction for a single 
qubit. If several qubits are coupled, the dephasing for 
the whole system is faster. Since the sources of dissipa- 
tion for different qubits are independent, the dephasing 
time gets N times shorter. In addition, the energy split- 
tings of the qubits should fit into the range provided by 
Eq. (||), and the levels within this range should be suffi- 
ciently different from each other to minimize the errors 
introduced by non-zero inter-qubit coupling between 2- 
bit operations (cf. Section ||) . This provides a restriction 
on the number of qubits and the number of operations 



TABLE I. Examples of suitable sets of parameters. 



c 


Cj 


Ej 


T 


L 


1-bit 


2-bit 


(aF) 


(aF) 


(mK) 


(mK) 


( M H) 


oper-s 


oper-s 


2.5 


400 


100 


50 


3 


8- 10 5 


650 


2.5 


400 


250 


125 


1 


8- 10 5 


500 


2.5 


400 


250 


125 


3 


8- 10 5 


1600 


40 


400 


40 


20 


1 


4- 10 3 


85 


10 


400 


100 


50 


1 


5- 10 4 


200 


40 


400 


100 


50 


0.5 


4- 10 3 


100 



The quantum measurement process introduces addi- 
tional constraints on the parameters. In order to demon- 
strate that the conditions assumed in this paper are real- 
istic we chose the charging energies -Eset, E qo and E- lnt 
as follows: The capacitance of the Josephson junction 
is Cj = 4.0 • 10~ 16 F, the gate capacitance of the qubit 
C = 2.5 • 10~ 18 F, the capacitances of the normal tunnel 
junctions of the SET C = 2.0 • 10~ 17 F, the gate capaci- 
tance of the SET C g = 2.5 • 10~ 18 F, and the capacitance 
between the SET and the qubit C int = 2.5 • 10" 18 F. 
We obtain: E SET ~ 20K, E qh ~ 2.5K, E int w 0.25K. 
Taking n qh = 0.35, A SET = 0.15 and eV tI = 48K we 
get AE » 3K, E C b = £ S et(1 ~ 2A SET ) ~ 14K, and 
Vtr/2 — Ecb ~ 10K (for definitions see subsection IV D| ). 
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We also assume 2-rra — 0.1. The measurement time in 
this regime is t ms w 0.25 • 10 4 ?i/(fcB IK) w 1.8 • 10~ 8 s. 
For this choice of parameters we calculate t m i x using 
Eq. (|62j). Assuming first Ej = 0.1K we obtain t m i X w 
1.4- Wh/(k B IK) w 1.0 • 10" 6 s. Thus t mix /t ms ps 55 and 
the separation of peaks should occur much earlier than 
the transitions happen. Indeed, the numerical simulation 
of the system ([l8| ) for those parameters given above shows 
almost ideal separation of peaks (see Fig. ||). On the 
other hand, for Ej = 0.25K, and we obtain i m i x /i ms ~ 9. 
This is a marginal situation. The numerical simulation 
in this case (see Fig. ^ shows that the peaks, first, start 
to separate, but, later, the valley between the peaks is 
filled due to the mixing transitions. 



These numbers demonstrate that the quantum manip- 
ulations of Josephson junction qubits, as discussed in this 
paper, can be tested experimentally using the currently 
available lithographic and cryogenic techniques. We have 
further demonstrated that the current through a single- 
electron transistor can serve as a measurement of the 
quantum state of the qubit, in the sense that in the case 
of a superposition of two eigenstates it gives one or the 
other result with the appropriate probabilities. 



P(m,t) 

0.005 




FIG. 9. P(m,t), the probability that m electrons have 
tunneled during time t. The parameters are those given in the 
text, Ej — 0.25K. The time is measured in nanoseconds. The 
initial amplitudes of the qubit's states: a — V0.75, b = \/0.25. 



P(m,t) 




FIG. 8. P(m,t), the probability that m electrons have 
tunneled during time t. The parameters are those given in the 
text, Ej — 0.1K. The time is measured in nanoseconds. The 
initial amplitudes of the qubit's states: a — V0.75, b — V0.25. 



B. Comparison with existing experiments 

The demonstration that the qubit is in a superposi- 
tion of eigenstates should be distinguished from another 
question, namely whether it is possible to verify that an 
eigenstate of a qubit is actually a superposition of two dif- 
ferent charge states, which depends on the mixing angle 
r\ as described by Eq. (|34|). This question has been ad- 
dressed in the experiments of Ref. Jll[ . They used a setup 
similar to the one shown in Fig. |4|Ja single-Cooper-pair 
box coupled to a single-electron transistor. They could 
demonstrate that the expectation value of the charge in 
the box varies continuously as a function of the applied 
gate voltage as follows from (|34|). Another experiment 
along the same lines has been performed by Nakamura 
et al. pjjj . They demonstrated by spectroscopy that the 
energy of the ground state and of the first excited state 
have the expected gate- voltage dependence of superposi- 
tions of charge states. 

Our theory can also describe stationary-state measure- 
ments of Ref. [ pT[ . The measurement was performed in 
the whole range of voltages including the vicinity of the 
degeneracy point, and we need a different pertrubation 
expansion which is valid in this region |2(i[ | . For this pur- 
pose we have to analyze the rates in the master equation 
([IT]) for general values of the mixing angle rj, relaxing the 
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requirement rj <c 1. To this end we rewrite the master 
equation (|49| ) in the qubit's eigenbasis (§) and consider 
the mixing terms as a perturbation. Then, the expansion 
parameter is EjEi nt / AE 2 (ri) which extends the region of 
the validity of the perturbative treatment. It turns out 
that each eigenstate of the qubit, |0) or |1), corresponds 
to a single, though ^-dependent tunneling rate To/i^). 
Thus, if the qubit is prepared in one of its eigenstates, 
only one peak is observed. This is the case even close 
to the degeneracy point rj ss ir/2 where the eigenstates 
are superpositions of two charge states with substantial 
weight of both components. On the other hand a charge 
state, being a superposition of two eigenstates, would 
produce two peaks in the current distribution. In this 
sense the measurement provides information about occu- 
pation probabilities of the eigenstates rather than charge 
states of the qubit. The correction to the tunneling rate 
^o/i (v) ~ r oc ± cos 2 rj is proportional to the average 
charge in the corresponding eigenstate of the qubit, in 
accordance with the experiments of Bouchiat et al. JlT| . 
Close to the degeneracy it is more difficult (takes longer) 
to distinguish the eigenstates since the peaks get closer. 

In our solution of the master equation we neglected 
dissipative effects due to the environment discussed in 
Section III. It is justified on time scales shorter than 



the environment-induced relaxation and dephasing times 
given by (p2f), (p3[). It is also justified at longer times as 
long as the environment-induced effects are weaker than 
the measurement- induced mixing, i.e. the rates (|2^), ( |23| ) 
are smaller than the mixing rate ( |62| ) . In this limit ther- 
mal relaxation is ineffective, and in the stationary regime 
(t — > oo) the qubit is in the equally-weighted mixture of 
two states, corresponding to the infinite effective temper- 
ature. The current in the SET, el 1 , is insensitive to the 
gate voltage in the qubit's circuit. In the opposite limit 
of weaker mixing, which is relevant to the experiments 
of Bouchiat et al. mixing is ineffective and thermal 
relaxation takes over. At low temperature k^T <C Ej 
dissipation keeps the system in the ground state, and the 
stationary current value is j-To(rj). 



C. Related theories 



It is also interesting to compare our proposal with the 
"quantum jumps" technique employed in quantum op- 
tics in general, and with the realizations of the qubits 
by trapped ions in particular (for a review see Ref. [Eq] ) . 
Indeed, the concepts are very close in spirit: the state 
of the system is examined by an external noncquilibrium 
current (electrons in our case and photons in the quan- 
tum jumps technique). There is, however, an important 
difference. In the quantum jumps measurements only 
one of the logical states scatters photons. Therefore, the 
efficiency of the measurement is limited by the ability to 
detect photons. In principle we could realize this situ- 
ation also in our system, if we bias the SET transistor 



such that different states of the qubit switch the transis- 
tor between the off and on regimes. Then the efficiency 
of the measurement is determined by the ability to detect 
individual electron — which is possible in single-electron 
devices, for instance by charging a single-electron box — 
and the measurement time would be given by the time it 
takes the first electron to tunnel. However, this mode of 
operation would require that the SET transistor is kept 
near the switching point, where thermal fluctuations and 
higher order processes could modify the picture substan- 
tially. Therefore, we have concentrated here on a situa- 
tion in which the SET transistor conducts for both states 
of the qubit, and the measurement requires distinguish- 
ing large numbers of charges or macroscopic currents. 
Accordingly, the measurement time is limited by the shot 
noise. 



D. Summary 

To conclude, the fabrication of Josephson junction 
qubits is possible with current technology. In these 
systems fundamental features of macroscopic quantum- 
mechanical systems can be explored. More elaborate de- 
signs as well as further progress of nano-technology, will 
provide longer coherence times and allow scaling to larger 
numbers of qubits. The application of Josephson junc- 
tion systems as elements of a quantum computer, i.e. with 
a very large number of manipulations and large number 
of qubits, will remain a challenging issue, demanding in 
addition to the perfect control of time-dependent gate 
voltages a still longer phase coherence time. We stress, 
however, that many aspects of quantum information pro- 
cessing can initially be tested on simple circuits as pro- 
posed here. 

We have further shown that a single-electron transistor 
capacitively coupled to a qubit may serve as a quantum 
measuring device in an accessible range of parameters. 
We have described the process of measurement by deriv- 
ing the time evolution of the reduced density matrix. We 
found that the dephasing time is shorter than the mea- 
surement time, and we have estimated the mixing time, 
i. e. the time scale on which the transitions induced by 
the measurement occur. 
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